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Abstract. - A detailed numerical and experimental study of the ordering of the low tem- 
perature tetragonal phase of CuAu is presented. The numerical simulations are based on a 
coarse-grained free energy derived from electronic structure calculations of CuAu, while the 
experimental results are obtained from in situ x-ray scattering. Both theoretical and experi- 
mental work indicate a subtle kinetic competition between the ordered tetragonal phase and 
the metastable modulated phase. 



What determines the microstructure of a typical alloy? While there is no simple answer to 
this question it is understood 0, |^ that a variety of morphological instabilities can be encoun- 
tered in non-equilibrium processing which create spatial inhomogeneities that are magnified 
during coarsening. In this paper experimental and theoretical techniques are combined to 
investigate the intricate morphologies that emerge during the sublattice ordering of the low 
temperature tetragonal phase of CuAu. The results of this study indicate a subtle competition 
between a metastable modulated phase and the stable low temperature ordered phase. At 
shallow quench temperatures or early times the modulated state is favored while at lower 
temperature or later times the ordered phase dominates. The results are supported by in situ 
x-ray scattering experiments and numerical simulations. 

The phases of CuAu that are relevant to this study are a high temperature disordered 
face-centered cubic phase, an intermediate temperature modulated orthorhombic phase and 
a low temperature ordered tetragonal phase. The modulated phase consists of a periodic 
array of antiphase domain walls arranged along one of the directions perpendicular to the 
tetragonal direction. The wavelength of the modulated phase is 10 underlying fee unit cells. 
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The transitions from the ordered phase to the modulated phase and from the modulated 
phase to the disordered phase are both first order and occur at Tqm = 658if and Tmd — 
683i4r respectively. Rapid temperature quenches through these transitions can lead to many 
interesting morphologies. In this paper quenches to the low temperature ordered phase (i.e., 
T < Tom) are considered. Earlier experiments]^ suggested homogeneous nucleation for 
quenches just below Tqm and spinodal ordering for deep quenches. 

In principle the ordering kinetics involve both structural relaxation of the lattice and a 
local rearrangement of atoms into the appropriate sublattice. Although it is desirable to use a 
microscopic model that describes these two distinct processes |j, it is difficult to explore the 
length and time scales of interest in this paper with such an atomistic approach. To study these 
long length and time scales a mean-field continuum model|^, |^ will be used. The mean-field 
model is based on a microscopic Hamiltonian derived from effective medium theory(EMT)|Q|. 
EMT is a semi-empirical approach based on the concept of an electron-density dependent 
contribution to the total energy of a system of atoms. The parameters of this model were fit 
to experimental values of the lattice parameter, bulk modulus and shear modulus of Cu and 
Au and to the heat of formation of the 50-50 alloy Q. The Hamiltonian describing Cu-Au 
alloys through the whole concentration range is defined by the EMT functional form and these 
parameters. The mean-field model is obtained from this Hamiltonian and describes the free 
energy of the system as a functional of the sublattice concentration and explicitly assumes that 
the structural relaxation is instantaneous. The model predicts the existence of the phases and 
transitions described above. As with most mean-field theories the predictions for the transition 
temperatures (i.e., Tmd = 1263K and Tqm — 1211iir) are not in good agreement with the 
experimental values |Q. An enhanced mean- field theory could be constructed to improve the 
quantitative details, but the competition between ordered and modulated phases predicted by 
the current model already agrees qualitatively with experiment. 

Assuming that the kinetics are driven by minimization of this free energy the following 
kinetic model can be proposed, 

drj/dt = -2T{27r/af{[aT + e(Vi - Vf,) + f\/i]rj - 4u7j^ + 6vtj^) + (, (1) 

where, {Cif,t)C{f,t)) = -2kBTr5{f~ f')S{t - t'), ar = 0.042(T - l2l9K)meV/K and r/ rep- 
resents the sublattice concentration. All parameters, except F, the mobility, were determined 
by the EMT calculations described above and are given in previous work[^. Rather than 
introducing a phenomenological or empirical form for T all times will be recorded in terms of 
a Ginzburg-Landau time defined such that the fastest growing Fourier mode grows as e*/*"^^ . 
For this model toL — l/[2r(27r/a)''(e^/4/ — ar)]- The subscripts _L and || refer to directions 
perpendicular and parallel to the ordering directions respectively. The interesting structures 
develop in the plane perpendicular to the ordering direction and for computational efficiency 
gradients in the 1 1 direction will be ignored. This model diff'ers from the standard kinetic model 
of sublattice ordering (Model A) since the coefficient of the V^^ term is negative. Eq. (||) was 
simulated on a discrete lattice using Euler's method to evaluate the time derivatives and a 
'spherical Laplacian' introduced by Oono and Puri[H to account for the spatial gradients. 
Relatively small space and time steps and large system sizes were used to eliminate spurious 
numerical effects. Eq. (|^) was simulated on system of dimension 2048 x 2048 grid points 
corresponding to a physical dimension of 0.5fim x O.bfim. 

The three equilibrium phases, disordered, modulated and ordered are described by (ry) = 0, 
r] ^ A sui(qof ■ n) and {rj) = rjo respectively, where n is an arbitrary unit vector in the 
ordering plane and Qo is the modulation wavevector. This differs from CuAu for which 
the modulation wavevector appears only along the Bravais lattice directions in the ordering 
plane. Transmission electron microscope workp, nol has indeed shown that the morphologies 
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in the modulated regime are quite anistropic. Nevertheless, even at this level of simplification 
the model contains many complexities that are important to the ordering process. Most 
importantly it incorporates the subtle competition between the modulated, disordered and 
ordered phases. 

To examine the kinetics of ordering below Tqm several simulated quenches were examined. 
In the first set of quenches the system was equilibrated in the disordered phase (i.e., at 
T — Tmd = 20-?^) and then instantaneously quenched to temperatures below Tqm- A 
qualitative picture of the subsequent dynamics is shown in fig. (|^) for quench temperatures of 
Tom — T = 5K and 15K . At both temperatures the very early stages are dominated by highly 
interconnected structures that contain many defects. The subsequent kinetics are dominated 
by a coarsening in which the average length scale of the patterns increases. This coarsening 
occurs mainly near defects in the pattern. Figure (|^) also indicates that the lower quenches 
coarsen more rapidly. This is due to the increase in the free energy difference between the 
modulated and ordered states at lower temperatures. Although not shown here this effect was 
verified at other quench temperatures (i.e., T — Tqm = 13K and 65K). 

To further illustrate the kinetics of ordering and for comparison with the experimental x-ray 
scattering patterns it is useful to consider the spherically averaged dynamic structure factor, 
S{q,t). S{q,t) is defined as, 

siq,t)^ J2 E 1' (2) 

where ri{q, t) is the discrete Fourier transform of ri{x, t). The ordered and modulated phases are 
represented by a peak(s) in S{q,t) centered around qo — and (Jm = ±.18/ A respectively. To 
examine the competition between the two phases the structure factor was fit to two symmetric 
peaks: an ordered peak, So{k,t), and a modulated peak, SM{k,t), centered around qo and 
qM respectively. For comparison with experiment the integrated intensities (i.e., J dqSoiq,t) 
and J dqSM{q,t)) were evaluated for each peak. Ideally, these quantities reflect the volume 
fraction of the two phases. In practice there is no clear distinction between 'ordered' and 
'modulated' regimes as is illustrated in morphological patterns shown in fig. (jl|). Nevertheless 
these quantities are useful for comparison with experiment as they are relatively insensitive to 
experimental resolution (which is not the case for peak heights and widths). The dynamics of 
the integrated intensities and some representative structure factors are shown in fig. (^. 

The structure factors in fig. (||) indicate that satellite peaks are generated at early times 
for quenches just below Tom consistent with the highly interconnected patterns shown in fig. 
(0a) and (|^c). The metastability of the modulated phase is highlighted by the persistence of 
the satellite peaks at late times as seen in fig. (j|a) and (||b). At the lower quench temperature 
(see figs. (||c:) and (^)) the central peak associated with the ordered phase dominates more 
rapidly. 

To test the numerical predictions and to provide additional insight into the ordering, x-ray 
scattering studies on the modulated and ordered superlattice peaks near (110) were carried out. 
The experiments were conducted on beamline X20C at the National Synchrotron Light Source. 
Several kinds of CuAu samples were studied - including bulk single crystals and polycrystalline 
films approximately 10/im thick. The detailed ordering kinetics differed slightly between 
samples but the qualitative features have been reproduced consistently. The results presented 
in this paper are from the polycrystalline films. The samples were held in a He atmosphere 
and in situ x-ray diffraction patterns were recorded by a linear position-sensitive photodiode 
array detector every 3s. Quench rates were typically of the order 5.5 degrees/s. Similar to the 
numerical simulation the prequench samples were equilibrated at T — Tmd = SOX and then 
rapidly quenched to temperatures below Tqm- All x-ray intensities, I{q,t), were normalized 
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with the incident intensity and a hncar function was subtracted from each pattern to remove 
the background. 

The x-ray scattering results for quenches to Tqm — T — 5K and 15K are summarized in 
fig. (|^) and can be compared directly with the numerical results shown in fig. (^). In this 
comparison it should be noted that the numerical results are for correlations in the sublattice 
concentration (77) which, for example, peak at q = for the ordered phase (i.e., 77 = rjo). In the 
experiment the same ordered phase leads to a superlattice peak or Bragg reflection centered 
around 2'K/d where d is the (110) plane spacing. Both experiment and simulation show that 
satellite peaks are generated at early times and that the persistence of these peaks is related 
to the quench temperature. Although not shown here, it was also found, in both experiment 
and simulation, that the early stage kinetics are very similar above and below Tqm- 

One feature that was observed in the experiment and not in the simulations was the shift 
in position with time of the x-ray peaks. This motion arises from changes in both the 
lattice spacing and modulation wavelength. The fundamental (200) peaks (not presented 
here) directly show the evolution of the lattice structure |T^. Examination of these indicates 
that, at the temperatures and time scales discussed below, the lattice has already relaxed to 
its tetragonal shape. Another significant effect is the shifting of the satellite peaks away from 
the central superlattice peak. These changes are typically quite large and occur in the early 
to intermediate time scales. For example in the quench to Tqm — T = IbK the modulation 
wavelength decreased from approximately 15.6 to 11.5 lattice constants on a time scale of 300s. 

Diffuse scattering studies of the disordered alloy have found satellites around the superlattice 
peaks with a "modulation" wavelength approximately 40% larger than in the equilibrium 
modulated state|l2j. In the early stages of ordering following a quench, these satellite peaks 
grow rapidly in height before moving outward to the equilibrium positions. 

The ordering kinetics can be strongly altered by changing the prequench state from a 
disordered phase to a modulated phase since, in the model, the modulated phase is metastable 
at low temperatures and the disordered phase is unstable. To examine this phenomenon two 
deep quenches were considered, one from a disordered initial condition (as before) and one 
from a modulated state. The numerical prequench state was created by seeding the system 
with large modulated regimes and evolving the system until no disordered patches remained. 
The resulting pattern was a modulated state containing many defects. For comparison the 
prequench experimental system was equilibrated in the modulated phase at T — Tqm = Sif . 

Quenching these states to Tqm — T = &hK led to the results illustrated in fig. (|^). The 
influence of the prequench state can be clearly seen in figs, (gi) and (|b) which compare 
configurations at the same time following the quench. The persistence of the prequench 
modulation is quite apparent in fig. (^). The erosion of this modulated phase occurs only at 
defects in the pattern. In principle the ordered regions can spontaneously nucleate in a purely 
modulated regime, however it appears that the nucleation rates are very small. For example, in 
one test simulation conducted at Tqm — T = Q^K no nucleation events were observed for times 
upto t/tcL = 5,000. This is also apparent although less obvious in the quenches described 
earlier (see fig. (0)). 

The experimental and theoretical results for these quenches are compared in figs. (^) and 
(Ijd). In fig. (^) the scattering intensity is shown at 1650 seconds after the quench for both 
prequench states. A similar comparison is provided in fig. (^) for the numerical simulations. 
In both experiment and simulation growth of the ordered state is seen to be strongly reduced 
for the modulated prequench state. This can be seen by comparing the amplitudes of the 
central and satellite peaks for the two different prequench states. 

In summary, the combined experimental and theoretical work presented in this paper has 
provided a detailed description of the kinetics of ordering in the low temperature tetragonal 
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phase of CuAu. This description highUghts the influence of the metastable modulated phase on 
the ordering process and the interesting kinetics that result. The continuum model provided 
a good qualitative prediction of the competition between the ordered and modulated phases 
as a function of time and pre- and post-quench temperatures. 
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Fig. 1. - Transient morphologies for Tqm — T = 5K and 15K obtained from numerical simulations. 
The light and dark regions respectively correspond to 77 < and 77 > 0. Only one twenty fifth of the 
full simulations cell (0.5/im x 0.5/im) is shown here. Figures (a) and (b) correspond to Tqm —T — 5K 
at times t/toL = 112 and 375 respectively. Figures (c) and (d) depict r; for Tqm — T= 15K at times 
t/toL =49 and 146 respectively. 



Fig. 2. - In figures (a) and (c) the rmrricrical simulated integrated intensities arc shown as a function 
of time for quenches to Tqm — T = 5K and Tqm — T= 15K respectively. The solid, long-dashed and 
short-dashed lines correspond to the modulated, ordered and total integrated intensities respectively. 
In figures (b) and (d) sample structure factors are shown at Tqm — T = 5K and Tqm — T = 15K 
respectively at the times indicated by the arrows in (a) and (c). The lines from bottom to top at 
q = correspond to the earliest to latest times. 



Fig. 3. - In figures (a) and (c) the experimentally measured integrated intensities are shown as a 
function of time for quenches to Tqm — T = 5K and Tqm —T= 15if respectively. The solid, long- 
dashed and short-dashed lines correspond to the modulated, ordered and total integrated intensities 
respectively. In figures (b) and (d) sample scattering patterns are shown at Tqm — T = 5K and 
Tom —T = 15K respectively at the times indicated by the arrows in (a) and (c). 



Fig. 4. - In figs, (a) and (b) configurations are shown for at t/toL ~ 18.8 for quenches from disordered 
and partially modulated phases respectively. As with fig. (1) only one twenty fifth of the full simulation 
cell is shown. In fig. (c) the experimental scattering intensity is shown at 1650s after the quench for 
both disordered (dashed line) and modulated (solid) initial states. A similar comparison is shown for 
the numerical structure factor in (d) at t/tcL = 18.8 after the quench. 
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